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Abstract — We present an algorithm to construct quantum 
circuits for encoding and inverse encoding of quantum convo- 
lutional codes. We show that any quantum convolutional code 
contains a subcode of finite index which has a non-catastrophic 
encoding circuit. Our work generalizes the conditions for non- 
catastrophic encoders derived in a paper by Ollivier and Tillich 
I quant-ph/0401134 ( which are applicable only for a restricted 
class of quantum convolutional codes. We also show that the 
encoders and their inverses constructed by our method naturally 
can be applied online, i. e., qubits can be sent and received with 
constant delay. 

I. Introduction 

Similar to the classical case a quantum convolutional code 
encodes an incoming stream of quantum information into an 
outgoing stream. A theory of quantum convolutional codes 
based on infinite stabilizer matrices has been developed re- 
cently, see [12]. While some constructions of quantum con- 
volutional codes are known, see [2], [3], [1], [5], [6], [7], 
[11], [12], some very basic questions about the structure of 
quantum convolutional codes and their encoding circuits have 
not been addressed so far, respectively have been addressed 
only in special cases. In this paper we focus on the question 
of which quantum convolutional codes have non-catastrophic 
encoders, respectively inverse encoders. 

Recall, that classically a code encoded by a catastrophic 
encoder has the unwanted property that — after code word 
estimation — a finite number of error locations can be mapped 
by the inverse encoder to an infinite number of error locations. 
For classical convolutional codes it is well-known that the non- 
catastrophicity condition is a property of the encoder and not 
of the code itself. Indeed, every convolutional code has both 
catastrophic and non-catastrophic encoders and therefore the 
choice of a good encoder is very important. 

In this paper we address the analogous question whether any 
quantum convolutional stabilizer code has non-catastrophic 
encoders and encoder inverses. Here the condition to be non- 
catastrophic has been shown in [12] to be that it has a constant 
depth encoder whose elementary quantum gates can be ar- 
ranged in form of a "pearl necklace", i. e., a regular structure in 
which blocks are only allowed to overlap with their neighbors 
with possibly some blocks spaced out. Furthermore, in [12] 
some conditions on the code have been given under which 
a non-catastrophic encoder exists. However, these conditions 



are quite strict and not applicable to an arbitrary quantum 
convolutional code. 

Using the matrix description of quantum convolutional 
stabilizer codes and transformations on this matrix which 
preserve the symplectic orthogonality, we show that a normal 
form can be achieved which corresponds to a very simple 
convolutional code. Reducing the dimension of this code 
by only a bounded factor, we obtain an even simpler code 
allowing online encoding and decoding. Furthermore, from 
the sequence of transformations one can read off a non- 
catastrophic encoder for a subcode of the original code whose 
dimension is reduced by the same factor. Asymptotically, the 
rate of the subcode and the original code are the same. 

II. Quantum Convolutional Codes 

Quantum convolutional codes are defined as infinite versions 
of quantum stabilizer codes. We briefly recall the necessary 
definitions and the polynomial formalism to describe quantum 
convolutional codes which was introduced in [12]. 

Definition 1 (Infinite Pauli Group): Let 

be the (real version of the) 2x2 Pauli matrices. Consider an 
infinite set of qubits labeled by the nonnegative integers N. 
Let M e {X, Y, Z} be a Pauli matrix. We denote by M; the 
semi-infinite tensor product /2 ® . . . ® /2 ®M(g)/2 (8) . . ., where 
AI operates on qubit i and I2 denotes the identity matrix of 
size 2x2. The group generated by all and Zi for i € N 
is called the infinite Pauli group Voo- For an element A — 
Ai (g) A2 <E) ■ ■ ■ S Voo the positions in which Ai is not equal 
to ±/2 is called the support of A. 

In the theory of block stabilizer codes, the elements of the 
Pauli group are labeled by tuples of binary vectors. Similarly, 
we can label the elements of the infinite Pauli group by a tuple 
of binary sequences, each of which is represented by a formal 
power series. Hence we get the correspondence 

i^iyX^Z^ := (-l)^(g)X"*Z^' 

\£>0 l>0 I 



where c G F2 and a = J2e>a '^tD^ and /3 — '}Zi>o l^tD^ ^re 
formal power series with coefficients in F2. In this representa- 
tion, multiplication of elements of Voc corresponds to addition 
of the power series. Furthermore, shifting an element A G Voo 
one qubit to the right corresponds to the multiplication of the 
power series by D. As we also allow to shift the operators 
by a bounded number of qubits to the left, we use Laurent 
series instead of power series to represent the elements of 
Voo- An element A G Voo with finite support corresponds to a 
tuple of Laurent polynomials. Recall that the field of Laurent 
series in the variable D with coefficients in F2 is denoted by 
F2 ((£»)) and recall further that it contains the ring F2[D, D"^] 
of Laurent polynomials. 

We are interested in shift invariant abelian subgroups of Voo, 
more specifically in those subgroups which can be generated 
by a finite number of elements and their shifted versions. 
The following definition introduces a shorthand notation for 
describing such subgroups. 

Definition 2 (Stabilizer Matrix): Let S be an abelian sub- 
group of 7^00 which has trivial intersection with the cen- 
ter of Vqo- Furthermore, let {gi, (72, ■ • • , <?r} where gi = 
{-lY'Xo^^Zp, with c, G {0,1} and K,/3,) G F2((D))" x 
^2{{D))^ be a minimal set of generators for S. Then a 
stabilizer matrix of the corresponding quantum convolutional 
(stabilizer) code C is a generator matrix of the (classical) 
additive convolutional code C C F2((i:'))" x F2((i:>))" 
generated by {ai,(3i). We will write this matrix in the form 



S{D) ^ {X{D)\Z[D)) 



1 Oil 


/3l \ 




fir I 



GF2((i?)) 



rx2n 



(1) 

In what follows we are only interested in those stabilizers 
which have a finite description. Hence we will consider only 
such stabilizer matrices {ij in which all entries are actually 
rational functions, i.e., elements of ¥2(0) . Eventually, we 
will require that all entries have finite support and are hence 
polynomials. 

Alternatively to Q a quantum convolutional code can also 
be described in terms of a semi-infinite stabilizer matrix S 
which has entries in F2 x F; 
matrix is as follows: 

/ Go Gi 
Go Gi 
Go 



S 



Gi 



The general structure of the 
\ 

(2) 
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G,i 







V : ■ ■ / 

The matrix 5* has a block band structure where each block is 
of size [n — k) X (m + \)n. All blocks have equal size and 
are comprised of m + 1 matrices Gq , Gi , . . . , Gm which are 
of size {n — k) X n each. In the second block, these to + 1 
matrices are shifted by n columns, hence any two consecutive 
blocks overlap in (m — l)n positions. 

Similar to the classical case the link between the polynomial 
description of eq. Q and the semi-infinite matrix eq. (|2|l is 
given by S{D) X]"=o GiD^. The band structure of eq. (|2j 



implies that for every qubit in the semi-infinite stream of 
qubits, there is a bounded number of generators of the stabi- 
lizer group that act non-trivially on that position. Moreover, as 
these generators of the stabilizer group have bounded support, 
their eigenvalues can be measured when the corresponding 
qubits have been received. Therefore, it is possible to compute 
the error syndrome for the quantum convolutional code online. 

Writing the stabiHzer in the form S{D) = {X{D)\Z{D)) as 
in eq. Q, it was shown in [12] that the condition of symplectic 
orthogonality of the semi-infinite matrix S can be expressed 
compactly in the form 



X{D)Z{1/DY + Z{D)X{llDf = 0. 



(3) 



On the other hand, we can start with an arbitrary self- 
orthogonal additive convolutional code over F2 ((!?))" x 
F2((£'))" to define a convolutional quantum code. In general, 
the generator matrix for such a code may contain rational 
functions, but there is always an equivalent description in 
terms of a matrix with polynomial entries [9]. The following 
theorem shows that for self-dual convolutional codes, all 
entries of a systematic generator matrix are in fact Laurent 
polynomials. 

Theorem 3: Let S{D) = {X{D)\Z{D)) with X{D) = I 
be a stabilizer matrix of a self-dual additive convolutional 
code over the rational function field V2{D). Then Z{l/D) = 
Z{Dy and all entries of Z{D) are Laurent polynomials. 

Proof: From condition (|3jl it follows that the code is 
self-dual if and only if Z{l/D) = Z{DY. Assume that 
Zij [D) is a proper rational function and not a Laurent poly- 
nomial. Then evaluating the series expansion of Zij{D) at 
1/D yields infinitely many negative powers. However, since 
Zji{D) contains only finitely many negative powers we get a 
contradiction. Hence all entries of Z{D) have to be Laurent 
polynomials. 

The symmetry Z{l/D) = Z{DY additionally implies that 
the diagonal terms Za (D) are Laurent polynomials of the form 



Z,,(Z?) =^c,(Z?-^ + i?0■ 



(4) 



III. Shift- INVARIANT Clifford Operations 

We are interested in quantum circuits which encode a convo- 
lutional quantum code. Recall that the controlled-not (CNOT) 
maps |a;}|y) \x)\x ®y) and that the controlled-Z (CSIGN) 
gates maps \x)\y) i— > (— l)^'^|x)|?/) (see [10]). We want that 
errors which happen during the encoding do not be spread out 
too far. A particularly bad example of spreading errors is given 
by the cascade CNOToo = UZo CN0T('''+^) where gates 
with smaller index i are applied first. The cascade CNOToo 
maps the finite support element X (g) /2 18) /2 ® • ■ • to the infinite 
support element X^X^X®... On the other hand the infinite 
cascade CSIGNoo = UZo CSIGN'^'^+i) does not have this 
behavior: indeed, a Pauli matrix Xi is mapped to Zi-iXiZi+i 
by this cascade and since it furthermore commutes with all Z 



operators, this shows that it maps finite support Pauli matrices 
to finite support Pauli matrices. The reason for this difference 
is that the sequence CSIGNqo can be parallelized to have 
finite depth (actually depth 2), whereas this is not possible for 
CNOToo- Clearly, any circuit of constant depth only leads to 
a local error expansion, i. e., Pauli matrices with finite support 
get mapped onto Pauli matrices with finite support. This gives 
rise to the following definition: 

Definition 4 (Non-catastrophic encoder): Let C be a quan- 
tum convolutional code and let £ be an encoding circuit for 
C. Then £ is called non-catastrophic if the gates in £ can be 
arranged into a circuit of finite depth. 

In the following, we consider infinite cascades of gates from 
the Clifford group that can be realized by quantum circuits 
with constant depth. Since the generators for the quantum 
convolutional code are obtained by shifting a fixed block an 
infinite number of times, we have to impose a shift invariance 
condition on any Clifford gate that we intend to apply to the 
code. This means that whenever a gate is applied it has to 
be applied also in a shifted version by an offset of n qubits. 
Similar to the approach in [8], the action of such operations 
on elements of the infinite Pauli group can be described as 
Unear transformations on the stabilizer matrix. As an example, 
the action of an infinitely replicated Hadamard gate H on a 
qubit is described in its action on the vectors {f{D),g{D)) G 

¥2(0)'^ by the matrix H = (^^^ ^ ^ since H^XH = Z 

and ZH = X. Similarly, all infinitely replicated versions 
of Clifford gates which only operate within a block and do 
not connect qubits between shifted blocks, correspond to the 
usual matrices in the symplectic group Sp2„(F2). 

More interesting are those operations which connect differ- 
ent blocks which have been shifted in time. An example is a 
CNOT gate which operates on a qubit i (control) and qubit 
j (target), where qubit j has been shifted by £ blocks. Recall 
that shifting by £ blocks corresponds to multiplying by D^. 
In this case we obtain that CNOT gate maps the stabilizer 
vector {xi,X2\zi, Z2) ^ {xi, X2 + xiD^\zi + Z2D~^ , Z2), i.e., 
X errors are propagated into the future and Z errors into the 
past. Note that by applying a sequence of CNOT gates we 
can actually map {xi, X2\zi, Z2) {xi,X2 + f{D)xi\zi + 
f{l/D)z2, Z2), where f{D) G F2[-D] is an arbitrary polyno- 
mial. A summary of the gates used is shown in Table H] It is 
important to note that all the operations shown in Table H] can 
be parallelized to have constant depth. 

IV. Computing an Encoding Circuit 

In the following we describe an algorithm which operate on 
the stabilizer matrix Q in order to produce a new stabilizer 
which is in a simpler form. We can act in two ways: (i) 
by applying row operations using an invertible matrix over 
¥2[D, D^^]. Apart from possible shifts, this does not change 
the stabilizer group, i. e., up to a possibly new initial qubit 
sequence (of bounded length) the quantum code is unchanged. 
We can also apply (ii) column operations given by an arbitrary 
element of the Clifford group shown in Table|l] Before we state 



TABLE I 

Action of various Clifford operations. 



unitaiy gate I 


J 


matrix U 


\U cxp(27r/2j 


] e C2x2 p = 
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CNOT(^'-'+*"',j ^ j 
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0^ 





(mod n) CNOT = 




\0 


1 
D-^ I) 


CSIGN(^J+^"),i ^ j 




(1 
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D'\ 

D-^ 




(mod n) CSIGN = 



i^O 


1 
l) 


Pi := CSIGN('''+'^"' Jj^O Pi= \ 


1 j 



Conjugation of the stabilizer group cS by the unitary gate U corresponds to 
the action of the matrices U on the columns of the stabilizer matrix S{D) = 
{X{D)\Z{D)). 



the algorithm we recall the Smith normal form [9] of a matrix: 

Theorem 5: Let M{D) e FaiD]' "'" be an r x n polyno- 
mial matrix. Then there exist polynomial matrices A{D) G 
GLr(¥2[D]) and B{D) G GL„(F2[£']), both having determi- 
nant one, such that M{D) = A{D)r{D)B{D), where r(D) 
is the r X n matrix 

/ 71 (^) \ 

V IriD) ••• Oj 

where the diagonal elements (elementary divisors) ji G F2 [D] 
satisfy "filji+i for z = 1, . . . , r - 1. 

Note that the Smith form can be computed for any matrix 
over an Euclidean domain, including the ring ¥2[D, D^^] of 
Laurent polynomials (see, e.g., [4]). For this, define the degree 
of a Laurent polynomial / = Y^iLi c^D^ with c^g ^ ^ c^^ 

as\h~io\. 

We will also need an observation about matrices which have 
already been partially brought into Smith form and which 
contain Laurent polynomials as entries. 

Lemma 6: Let M{D) G ¥2[D, D-^Y"""- be a matrix con- 
taining Laurent polynomials and which has the form M = 
(diag{-f,{D))\U{D)), where U(D) G ¥2[D, D-^Y''^'^-''l 
Assume that for at least one i we have that 7^ does not divide 
the Laurent polynomials contained in the ith row of U{D). 
Then at least one of the polynomials 7j'(-D) arising in the 



Smith normal form of M (D) (after the denominators have 
been cleared by row-wise multiplication of powers of D) has 
a strictly smaller degree than the corresponding 'yi{D). 

Proof: Without loss of generality, we consider the first 
row (7i(D),0,...,0,/ip),...,/„_,p)) of Afp), where 
the fi{D) are Laurent polynomials. Clearing the denominators 
by a suitable power leaves us with D^^i{D) and polyno- 
mials fl{D) D^fi{D). Computing the Smith normal form 
we obtain the gcd of D'^ji{D), f[{D), . . . , f'^_^{D) which by 
assumption has to be a proper divisor of 71 (_D). ■ 
Next, observe that by using Clifford gates (acting on the 
X-part only) we can implement the matrix B{D) used in the 
Smith normal form. The reason for this is that in the computa- 
tion of the Smith normal form only elementary operations and 
permutations are necessary [9, Section 2.2]. We can realize 
these operations using the CNOT gates and permutations 
of the qubits, which can also be realized by CNOT. Left 
multiplication by an invertible matrix does not change the 
stabilizer, so there is no need to implement the matrix A{D) 
as quantum gates. 

Algorithm 7: Let a polynomial stabiUzer matrix S{D) — 
{X{D)\Z{D)) e Fail?]'''' 2" of full rank be given. 

1) Compute matrices A{D) and B{D) which realize the 
Smith normal form for X{D). Factor the matrix B{D) 
into elementary matrices of the form CNOT and per- 
mutations of qubits. Apply these operations to the code 
to obtain the new stabilizer matrix 

T{D) 




Zi{D) Z2{D) 



where T{D) is a diagonal matrix with non-zero polyno- 
mial entries of rank s and Zi{D) e ¥2[D, £)-i]''xs and 
Z2{D) G ¥2[D,D-^Y'^^'^-^^ are matrices with Laurent 
polynomials as entries. 
2) While the Z2{D) part of S{D) is not zero, repeat the 
following steps: 
• Use Hadamard gates H to swap Z2{D) into the X- 
part yielding 

T{D) 




S{D) 



X2{D) 



Zi{D) 



with X2{D) ^ Z2{D). 

• If y{D) has full rank and if all polynomials in 
row j of X2{D) are divisible by jj{D) for all 
j = 1, . . . , r, then use CNOT-gates to obtain zeros 
in both X2iD) and ^2(1?). 

• Else recompute the Smith normal form of the X- 
part and get either smaller elementary divisors or 
all polynomials in Z2{D) are multiples of the cor- 
responding elementary divisor. The degree of the 
elementary divisors decreases because of Lemma |6l 

3) The stabilizer matrix S{D) is now of the form S{D) = 
{T{D)0\Zi{D)0), where T{D) has a rational inverse 
since S{D) has full rank. 

4) From Theorem |3] it follows that all entries in the rows 
of Zi{D) are divisible (as Laurent polynomials) by the 



corresponding element of T{D) (consider the matrix 
T-^S{D) = (/ 0\T~^Zi{D) 0) which co ntains Laurent 
polynomials only). Hence, using CSIGN gates, clear all 
off-diagonal terms in Zi{D). 

5) From © it follows that we can cancel the diagonal of 
the matrix Zi{D) using the gates P and Pi. 

6) Finally, use Hadamard gates H to obtain Z-only gener- 
ators in diagonal form. 

This algorithm transforms the original stabilizer matrix into 
a stabilizer matiix Si{D) (0 0|r(L»)0) with r(L>) = 
diag(7i(Z))). In case ^i{D) — 1, the only possible sequence 
of states formed by the z* qubit of all blocks is |0)|0) . . . If 
^i{D) = D^, there are no constraints on the first £ qubits. 
Otherwise, the state |co)|ci)... corresponding to the power 
series expansion of l/ji{D) — J2e>o^i^^ ^^'^ shifted 
versions are allowed, too. As the sequence (q)^ is periodic, 
there are only finitely many different shifted versions. We 
ignore these additional states as they would require an infinite 
cascade of CNOT gates. As an example for this behavior 
consider the states |0)|0)... and allowed by the 

single qubit Z-generator 1+D. 

In case T{D) ^ I, which corresponds to catastrophic 
encoders in the classical case, we consider the code Co with 
stabilizer matrix Sq{D) := (0 0|/0). Now, Co is a proper 
convolutional subcode of the code Ci with stabilizer matrix 
Si{D). The dimension is only decreased by a bounded factor 
depending on r(D). In case T{D) = /, we have Cq — Ci. 

The subcode Cq has a very simple structure: a sequence of 
n — k qubits in the state |0) alternates with a sequence of k 
qubits Encoding for Co is done by inserting qubits in the 
state |0) into the input stream. To obtain a state of (a subcode 
of) the original convolutional quantum code C, apply the gates 
corresponding to the elementary matrices used in the algorithm 
in reversed order. The corresponding elementary gates are only 
Clifford gates which have to be replicated infinitely often. All 
elementary gates used can be parallelized into finite depth 
which implies that the operations can be carried out online. 
Hence we have shown the following result: 

Corollary 8: Let S{D) be the stabilizer matrix of a quan- 
tum convolutional code C. Then there exists a convolutional 
subcode Csub Q C with a non-catastrophic encoder and encoder 
inverse such that asymptotically the rates of Csub and C are 
equal. Moreover, the encoder and its inverse only use Clifford 
gates and allow for online encoding and inverse encoding. 

V. Example 

Consider the F4-linear rate- 1/3 convolutional code from [6, 
Table VI]) with generator matrix 



GiD) = {l + D I + ujD 1+ZuD). 
The corresponding stabilizer matrix is 

S{D)-- 



1 + D 1 1 + D 
ODD 
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1 + D 1 + D 1 
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Fig. 1. Encoding circuit for a rate 1/3 convolutional quantum codes. Every gate has to be repeatedly applied shifted by one block, i.e. three positions down. 
Note that the CSIGN gates are diagonal and hence can be arranged in any order. For the CNOT gates, each gate has to be repeated in its shifted version 
before the next gate can be applied. 



The first sequence of CNOT operations transforms the first 
row of X{D), and we obtain: 
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+D 
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(D^ + l)/D 
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Invertible row-operations do not change the stabilizer group, 
so adding times the first row to the second yields 
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+ D 
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D D 
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Again using CNOT, we transform the second row of X{D): 

(5) 
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Using CSIGN, we can clear the off-diagonal terms in the first 
row of Z{D), 
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and similar for the second row 
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Finally, using P and H we get Z-only generators: 
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To clear the entries in the first row of the Z-part of (|5jl, we 
need CSIGN gates whose target lie in the block before that 
of the target (see the third set of gates in Fig. 0. Therefore, 
encoding can only start in the second block. In the first block, 
all qubits are initialized to |0). Additionally, we ignore that the 
term D in the second row implies that there is no operator Zi 



in the stabilizer acting on the second qubit of the now second 
block and constrain the input to |0). 

A circuit for encoding is obtained by reversing the order 
of the transformations. The encoding circuit is illustrated in 
Fig. n Note that the circuit extends over three blocks, i.e., 
has total memory two. This is reflected by the fact that in this 
example the operations used by the algorithm to clear entries 
only involved Laurent polynomials of degree at most two. In 
contrast, an encoder for the classical convolutional code over 
F4 given by G{D) can be realized with total memory one. 

VI. Conclusion 

We have shown that the quantum convolutional codes 
obtained from self-orthogonal classical convolutional codes 
always have a subcode which has asymptotically the same 
rate and allows for non-catastrophic encoders. This shows that 
errors affecting these codes do not propagate in an unbounded 
fashion during the decoding process. For simplicity, we have 
presented the algorithm for qubit systems only, but as in [8], 
the technique applies to non-qubit systems as well. 
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